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ABSTRACT 


•  t- 


It  is  shown  that  those  interactions  of  an  elastic  body  8  with  an 
elastic  environment  E  manifested  by  local  surfacial  loadings  should  be 
modelled  by  boundary  operators  featuring  a  live  ^system  of  surfacial  forces, 
i . . e . ,  a  vector  field  s,  defined  over  the  boundary  (  38  of  8  and 
representing  the  surfacial  force  per  unit  area  exerted  by  E  over  8  ,  which 
depends  functionally  in  a  non-trivial  way  on  the  deformation  u  of  8» 


In  particular,  under  the  assumptions  that  the  constitutive  equation  of  8 
is  compatible  with  ellipticity ,  and  s  is  a  function  of  the  appropriate 
restrictions  to  3$  of  u  and  its  gradient,  it  is  also  shown  that  it  is 
reasonable  to  require  that  the  resulting  live-boundary  condition  of  traction 
be  genuine,  i.e. ,  a  normality  condition  prevail  insuring  the  primary 
consistency  of  the  pair  of  field  and  boundary  operators. 


In  addition  to  the  analytical  difficulties  to  be  expected^  (-cf.  [7]/' 
Chapter  2)‘,  when  normality  does  not  hold  the  mechanical  interpretation  of 
boundary  conditions  fails  to  be  unique.  That  this  is  indeed  the  case  is 
demonstrated  by  producing  an  explicit  example  in  linearized  elastostatics  of  a 
boundary  operator  which  can  be  interpreted  as  a  non-genuine  live-boundary 
condition  of  traction  or  as  a  dead-boundary  condition  of  frictionless  contact. 
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SIGNIFICANCE  AND  EXPLANATION 


In  the  real  world,  most  loadings  "follow"  the  deformation  or,  in  the 
terminology  used  in  this  paper,  are  live,  in  that  their  direction  and/or  their 
magnitude  may  vary  with  the  deformation  itself  (the  hydrostatic  pressure  on  a 
submerged  object  and  jet  loading  are  common  examples  of  live  loadings? 
buckling  of  elastic  rods  and  plates  and  panel  flutter  are  examples  of  statical 
and  dynamical  problems  from  structural  engineering  and  aeroelasticity , 
respectively,  where  live  loadings  make  an  essential  feature). 


In  spite  of  its  practical  importance  and  deep  mathematical  interest,  a 
general  treatment  of  live  loadings  in  linear  and  non-linear  (and,  especially, 
three-dimensional)  elasticity  is  still  lacking. 


The  purpose  of  this  paper  is  to  delineate  circumstances  under  which  one 


can  expect  a  given  boundary  condition  to  represent  realistically  an  unambigu¬ 
ous  mathematical  model  for  live  surface  loading.  It  is  shown  that  the 
ellipticity  up  to  the  boundary  of  the  field  operator,  as  well  as  the  normality 
of  the  pair  of  field  and  boundary  operators,  both  play  a  crucial  role 
(ellipticity  is  a  requirement  of  physical  and  mathematical  plausibility  of 


general  use  in  elasticity?  normality  is  a  basic  mutual  consistency  requirement 

for  field  and  boundary  operators  in  the  theory  of  elliptic  systems,  receiving 

’cession  For 

a  precise  mechanical  interpretation  in  this  paper).  'JTIS  GRA&I  ~ 
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THE  ROLE  OF  ELLIPTICITY  AND  NORMALITY  ASSUMPTIONS  IN 
FORMULATING  LIVE-BOUNDARY  CONDITIONS  IN  ELASTICITY 

P.  Podio-Guidugli* ,  G.  Vergara  Caffarelli**  and  E.  G.  Virga*** 

1 .  Introduction 

The  mechanical  interaction  of  a  body  with  its  environment  is  a  varied  and 
complex  phenomenon  which  is  difficult  both  to  observe  and  describe  carefully. 
Boundary  conditions  modelling  non-trivial  interactions  are  difficult  to 
formulate  mathematically?  besides,  they  may  easily  lead  to  unusual  boundary 
value  problems,  perhaps  even  so  much  unusual  as  to  challenge  and  defeat  an 
expert  analyst.  Not  only  that,  blunt  generalization  of  the  few  well- 
understood  types  of  boundary  conditions  may  pose  rather  subtle  problems  of 
interpretation  and  classification  to  the  expert  mechanist. 

In  Section  2,  the  bulk  of  this  paper,  we  develop  a  train  of  reasoning 
leading  to  formulate  an  unambiguous  boundary  condition  of  traction, 
reproducing  the  non-linear  interaction  of  an  elastic  body  and  an  elastic 
environment  in  strict  contact  with  it?  the  emphasis  here,  in  accordance  with 
our  title,  is  on  the  role  played  by  ellipticity  and  normality  hypotheses. 

In  Section  3,  we  briefly  discuss  the  importance  of  ellipticity  and 
normality  in  problems  of  elastostatics  linearized  about  a  stressed  equilibrium 
placement.  In  particular,  we  show  by  means  of  an  explict  example  that,  if  the 
environment  is  live  but  normality  fails  to  hold,  one  can  easily  run  into 
taxonomic  troubles. 
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2.  The  Genuine  Live-Boundary  Condition  of  Traction 

Let  a  (possibly  unbounded)  open  reg-ilar  region  8  of  ,  with  M  =  2 
or  3,  be  separated  by  its  exterior  E  by  a  smooth  orientable  surface  S , 
the  common  boundary  of  8  and  E ,  with  n  the  exterior  unit  normal  to  8 . 

Once  constitutive  prescriptions  have  been  given,  we  may  think  of  8  as 
the  body ,  and  of  E  as  the  environment  of  the  body  8 .*  Then,  given  any 
judiciously  chosen  set  of  data,  we  may  formulate  a  variety  of  mechanically 
significant  problems,  anyone  of  which  ultimately  reduces  to  determine  the 
state  ( i  .e .  ,  the  triplet  of  displacement,  strain  and  stress)  over  B  U  E, 
provided  a  suitable  set  of  jump  conditions  has  been  stipulated  for  the  state 
at  the  points  of  S. 

Loosely,  our  goal  is  to  sketch  a  credible  scenario  where  a  large  class  of 
such  problems  for  8  U  E  can  be  reduced  through  a  rational  procedure  to  one 
corresponding  problem  for  8  =  8  U  S ,  with  the  presence  of  E  replaced  by  a 
set  of  boundary  conditions  on  S.  it)  anticipate  one  of  our  main  findings, 
perhaps  not  about  to  come  as  any  great  surprise  to  the  knowledgeable  reader, 
these  boundary  conditions  will  as  a  rule  be  of  the  live  type,  i  .e . ,  they  will 
be  prescriptions  of  the  surfacial  stress  vector  which  depend  functionally  in  a 
non-trivial  way  on  the  solution  in  5. 

Now,  let  u,  T  denote  the  restrictions  to  S  of  the  unique  extensions 
u,  T  to  o  of  smooth  displacement  and  stress  fields  u,  T  over  o ; 
further,  let  (u,t),  with  t  =  Tn ,  denote  the  pair  of  fields  over  5  con¬ 
sisting  of  the  displacement  u  and  the  stress  vector  t;  finally,  let  (v,s) 


*  Notice  that  here  we  implicitly  choose  to  not  consider  the  interesting  case 
when  8  itself  is  thought  of  as  a  material  surface  (in  the  manner,  e .g. ,  of 
f'TJRTTN  s,  yuppoch  f  1  ]  ) ,  and,  therefore,  made  the  object  of  a  peculiar  constitu¬ 
tive  prescr i pt ion .  To  contemplate  such  a  case  would  require  a  more  delicate 
ana  1 ys i s . 
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be  defined  over  S  in  a  completely  analogous  way,  but  starting  from  displace¬ 
ment  and  stress  field  v  and  I  defined  over  E .  We  wish  to  interpret 
v(x),  s(x)  as,  respectively,  the  displacement  of  boundary  point  x  of  the 
environment  and  the  contact  force  exerted  at  x  by  the  environment  over  the 
body . 

A  rather  spontaneous  assumption  is  that  such  contact  force  depends 
functionally  on  the  displacement  field  over  E  U  S .  Suppose  that  the  stress 
response  of  E  to  deformation  have  a  local  character,  as  is  the  case  for 
an  elastic  environment.  Then,  for  x  e  S  fixed,  and  for  any  given  field  v 
on  E,  one  finds  natural  to  choose  the  functional  dependence  of  s  on  v  as 
follows 

(2.2)  s(x)  =  ^ ( v( x) , Dv(x ) ) , 
where  D  is  the  gradient  operator. 

A  slight  generalization  of  such  a  choice,  of  importance  in  our  further 
developments,  is  suggested  by  the  splitting  of  the  gradient  operator  D,  at  a 
point  x  eS,  into  its  normal  and  tangential  parts  Dn  and  Dt: 

(2.3)  Dv(x)  =  D  v(x)  +  D  v(x)  , 

n  r 

with 

(2.4)  D  v(x)  =  3  v(x)  0  n(x)  ,  D  v(x)  =  D  v(x)  , 

n  n  t  t 

and  with  3  denoting  the  operation  of  differentiating  along  the  direction 
n 

n(x),  the  normal  to  S  at  x.  Then,  in  view  of  (2),  (3)  and  (4),  one  is 
induced  to  replace  (2)  by 


(2.5) 


s ( x )  =  (5  (v(x)  ,Dfcv(x)  ,3nv(x) ) 


Our  next  step  is  to  lay  down  a  set  of  jump  conditions  for  the  state  at 
the  points  of  S . 

We  say  that  the  body  8  and  its  environment  E  have  a  strict  contact  in 
a  neighbourhood  of  a  point  x  of  their  surface  of  separation  S  if 

(2.6)  u(x)  =  v(x)  &  t(x)  =  s(x)  . 

In  a  strict  contact,  both  the  displacement  and  its  tangential  gradient 
suffer  no  jump  at  S  j  besides,  the  continuity  of  the  stress  vector  at  S 
implies  in  general  restrictions  on  the  jump  condition  for  the  normal 
derivative  of  displacement.  To  motivate  the  ad  hoc  assumptions  we  shall  make 
to  determine  those  restrictions,  we  turn  momentarily  to  linear  elasticity. 

Suppose  that  8  and  E  ,  besides  being  in  strict  contact,  are  made  up 
with  two  linearly  elastic  materials  with  elliptic  elasticity  tensors  C  and 
D,  respectively,  so  that  C,  say,  obeys 

( EC ) °  det  C°(a)  ?  0  for  c°ih*a*  “  Cijhkajak'  a  J*  0. 

If  C  (D)  can  be  firstly  extended  and  then  restricted  to  S  in  such  a  way  as 
to  preserve  ellipticity,  the  matrix  field 

(2.7)  A  =  C°(n)  (8  =  D° ( n) ) 

has  invertible  values  on  S.  Thus,  as 

(2.8)  t  =  (CDu)n  =  AS  u  +  (CD  u)n, 


we  have  that 


a  formula  which  shows  that,  once  u  is  assigned  on  S  ,  t  and  3nu  are 
essentially  interchangeable  bits  of  information  at  a  point  of  S.  Likewise, 

(2.10)  3  v  =  B~ 1 { s  -  (DD  u) n )  . 

n  t 

Thus,  in  particular,  if  8  and  E  are  in  strict  contact,  linearly  elastic 
and  elliptic,  the  normal  derivative  of  displacement  suffers  at  5  a  discon¬ 
tinuity  depending  on  the  local  values  of  the  stress  vector  and  the  tangential 
gradient  of  displacement,  as  (A,C)  and  (B,D)  may  be  expected,  in  general, 
to  differ. 

Remark  1  It  is  important  to  note  that  ellipticity,  and  not  linearity  of  the 
constitutive  law,  is  the  crucial  assumption.  Indeed,  for  B  comprised  of  a 
non-linearly  elastic  material,  (8)2  is  replaced  by 

(2.11)  t  -  T(3  u8n  +  D.  u)n  , 

n  t 

_  • 

and  this  equation  is  locally  solvable  for  3ru  if  the  Ellipticity  Condition 
prevails,  i .e. , 

3t 

(EG)  det  C(a)  t  0  for  Cih(a)  =  T.-i^a-ja^  ,  a  ^  0,  T  =  - —  . 

3  (Du) 

Remark  2  For  boundary- value  problems  ruled  by  PDE's,  ellipticity  hypotheses 
typically  involve  the  coefficients  of  the  field  operator;  here,  EC  plays  an 
important  role  at  the  boundary. 

In  the  light  of  the  above,  if  one  assumes  that  8  and  E  have  a  strict 
contact,  then  (6)  allows  one  to  give  (5)  the  following  provisional  form: 


(2.12) 


s(x)  *jj(u(x),  D^u(  x ) ,  S^vfx))  . 


It 


Moreover,  for  elliptic  elastic  pairs  of  a  body  8  and  an  environment  E  , 

•  • 

seems  safe  to  assume  that  both  3nu  and  3^v  depend  on  the  tangential 
gradient  and  the  stress  vector 

(2.13)  3  u  =  U (Du , t ) ,  3  v  =  v(D  u,s) 

n  t  n  t 

and  that  the  dependence  on  the  stress  vector  is  essential,  in  that  the 
3U  3U 

matrices  -r — /  t —  are  invertible. 

3t  3s 

By  (12)  and  (13),  there  exists  a  function  g  such  that 

(2.14)  s(x)  =  g(u{ x ) ,  Dtu(x),  s(x)J  . 


Last  step  towards  our  goal,  we  assume  the  condition  (14)  is  genuine,  in  that 
for  no  assignment  of  u,  Dfcu  it  does  reduce  to  an  identity  in  s. 
Accordingly,  we  require  that 


( QQ ) 


det  L  /  0,  with  L 


G£  implies  that  (14)  is  locally  solvable  for  sj  then,  s  can  be  expressed 
as  a  function  of  u  and  Dtu: 

(2.1 5)  s(x)  =4(u(x),  Dtu(x))  . 

The  last  formula  suggests  that  in  formulating  boundary  value  problems  the 


contact  force  exerted  by  the  environment  on  the  body  ought  chosen  in  general 


to  be  of  the  live  t\ 


Remark  3  The  special  dependence  (15)  of  s  on  boundary  information 
concerning  the  solution  has  been  considered  by  SPECTOR  [2,3],  who  has  termed 
simp.1  e  sue)]  sue  facial  loading  operators.  SPECTOR  has  motivated  his  choice  by 
the  observation,  due  to  GURTIN  [4],  that  pressure  loading,  a  well-known 
example  of  live  loading,  is  indeed  simple. 


Remark  4  For  two  linearly  elastic  elliptic  bodies  in  strict  contact,  (12) 
takes  the  form 

s  =  B3  v  +  (DD  u)n 
n  t 

and  is  readilly  inverted  as 

(2.10)  3  v  =  B~ 1 ( s  -  (DD  u)n )  , 

n  t 

which  in  turn  corresponds  to  (13)2  •  Therefore,  (14)  reduces  to  a  sheer 
identity,  and  Qg  cannot  be  satisfied.  This  indicates  that,  in  general,  it  is 
not  sound  to  require  GC  in  a  contact  situation.  However,  on  passing  from 
contact  to  boundary  conditions,  we  believe  that  one  can  regard  (13)2  as  a 
separate  constitutive  assumption,  independent  of  (12)  and  consistent  with  it 
in  the  peculiar  sense  made  explicit  by  GC. 

Substitution  of  (15)  into  (13)2  an<*'  when  use  is  made  of  (6)2,  into 
(13)^  allows  us  to  complete  with  a  suitable  jump  condition  for  the  normal 
derivative  of  displacement  the  set  of  jump  conditions  for  the  state  to  be 
stipulated  at  a  point  of  S.  In  summary,  under  hypotheses  of  strict  contact, 
ellipticity  and  genuinity,  the  displacement  u  and  its  tangential  gradient 
Dtu  must  be  continuous,  whereas  the  normal  derivative  S^u  of  displacement 
may  suffer  a  jump  described  by  a  function  of  u  and  Dtu,  whose  choice  has 
constitutive  character. 

We  are  now  in  a  position  to  state  a  boundary  condition  of  traction  which 
would  describe  a  large  class  of  contact  actions  exerted  by  an  environment  of 
elastic  type  on  an  elastic  body. 

Naively,  one  would  write  a  live-boundary  condition  of  traction  for  B  as 

t ( x )  =  s ( x )  for  X  e  s. 
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or  rather  (cf .  ( 2 )  and  (11)),  as 

(2.16)  X  (u,  Dfcu,  3ru)(x)  =  0  for  x  e  S, 

with  the  boundary  operator  X  defined  by 

®  •  A  © 

(2.17)  X(u,  Du,  3  u)  =  T(3  ufin  +  D.u)n  -  /(u,  3  u8n  +  D.u)  . 

t  n  n  t.  u  n  u 

On  the  other  hand,  in  view  of  (6),  (13)^  and  (15),  the  boundary  operator 
should  have  the  following  form 

(2.18)  *(u,  Du,  3  u)  =  3  u  -  a(u,  Du)  . 

t  n  n  t 

Now,  X  is  locally  reducible  to  £  if  the  implicit  equation 

X(u,  Du,  3  u)  =  0 
t  n 


_  • 

can  be  locally  solved  for  °nu;  this  is  the  case  if  the  following  Normality 
Condition  holds: 


(SC) 


det  M(n)  f  0  for  M(n)  =  C(n) 


fn,  f 


3  (Du) 


Within  the  scenario  we  have  set  up,  accepting  is  mandatory.  We 

a 

observe  that  EC  implies  NC  when  the  environment  is  dead,  and  ^  reduces  to  an 

A 

assigned  function  of  x  only  over  ;  more  generally,  when  fa  is  indepen¬ 
dent  of  the  deformation  gradient;  finally,  when  is  simple  in  the  sense  of 
Remark  3. 

We  call  (16),  with  X  restricted  by  NC,  the  genuine  live-boundary 


condition  of  traction. 


3 .  Ellipticity  and  Normality  in  Linearized  Elastostatics 

When  one  passes  from  a  local  study  of  a  non-linear  problem  to  the 
corresponding  linearized  problem,  both  Eg  and  tig  continue  to  play  an  important 
role.  Indeed,  our  present  tig  is  a  version  appropriate  to  elasticity  of  the 
homonymous  condition  introduced  by  GEYMONAT  [5]  and  GRUBB  [6]  in  the  general 
theory  of  linear  elliptic  systems,  and  appears  to  be  the  natural  extension  to 
systems  of  the  normality  condition  for  scalar  operators  (vid.  e.g.  [7],  chapter 
2). 

If  gg  and  gg  prevail,  we  have  shown  in  [8]  that  a  Green  formula,  which 
generalizes  the  reciprocity  formula  of  Betti  of  classical  elasticity,  can  be 
associated  to  the  traction  boundary  value  problem  of  linearized  elastostatics 
with  live  loads.  Such  a  Green  formula  allows  one:  (i)  to  define  formal 
ad joints  to  both  the  field  and  boundary  operator;  (ii)  to  determine  a  set  of 
conditions  sufficient  for  self-adjointness;  (iii)  to  state  compatibility 
conditions  on  the  data  necessary  and  sufficient  [9]  for  solvability  of  the 
underlying  boundary-value  problem  in  a  familiar  Hilbert  space  setting. 

We  now  wish  to  give  an  example  of  the  ambiguities  t.  at  may  accompany  the 
failure  of  Ng. 

In  view  of  (2.16),  a  linearized  live-boundary  condition  of  traction  can 
be  written  as 

(3.1)  (TDu( x) )n( x)  -  fDu(x)  -  Fu(x)  =  0, 

where  u,  with  slight  abuse  of  notation,  is  now  used  for  the  first  approxi¬ 
mation  of  the  displacement  from  the  reference  placement; 

(3.2)  t ( x )  =  ( TDu ( x ) ) n ( x ) ; 

(3.3)  s(x)  =  f Du ( x )  +  Fu( x ) ,  with  F  =  . 
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Choose  ft  to  describe  a  live  environment  such  that 
’  '  fihk  ~  pilTljhknj'  Fih  =  “pih  ' 

n  P  a  perpendicular  projection  field  over  S.  Then,  it  is  easy  to  check 
r  the  matrix  M(n),  which  appears  in  the  statement  of  reduces  to 

M(n)  =  (t  -  P)C(n), 

if.,  therefore,  singular,  whereas  (1)  takes  the  aspect 

b  >  (  1  -  P )  t  +  Pu  =  0 , 

e.psj  ’  alently  , 

7)  (1-P)t=0  &  Pu=0. 

promptly  recognizes  in  (7)  a  generalized  boundary  condition  of  the  dead 
>■  i:;‘.  ■  10]  ,  Chapter  V).  In  particular,  when  P  =  nBn,  (7)  reduces  to  the 
(  il’.i  !  contact  boundary  condition: 

t  -  (t*n)n  =0  &  u*n  =  0  . 

f,  I  as  in  (4)  and  P  =  n8n,  the  boundary  condition  (1)  can  be 
>  i  1  1  dasrified  as  a  live-boundary  condition  of  traction  or  as  a  dead- 

.  :  i:  ■..inlit.u'ni  of  contact. 

■O'*.;  similar  examples  can  be  concocted.  But,  the  validity  of  NC 

t  lu  i''-  any  pathology. 
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''Vymor.at. .  We  also  wish  to  thank  W.  Williams  and  M.  E.  Gurtin  for  many 
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